Mapping the gravitational wave background 



Neil J. Cornish| 

Department of Physics, Montana State University, Bozeman, MT 59717, USA 

Abstract. The gravitational wave sky is expected to have isolated bright sources 
superimposed on a diffuse gravitational wave background. The background radiation 
has two components: a confusion limited background from unresolved astrophysical 
sources; and a cosmological component formed during the birth of the universe. A map 
of the gravitational wave background can be made by sweeping a gravitational wave 
detector across the sky. The detector output is a complicated convolution of the sky 
luminosity distribution, the detector response function and the scan pattern. Here we 
study the general de-convolution problem, and show how LIGO (Laser Interferometric 
Gravitational Observatory) and LISA (Laser Interferometer Space Antenna) can be 
used to detect anisotropics in the gravitational wave background. 
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1. Introduction 

Gravitational wave detectors are fairly blunt instruments. In the low frequency limit, 
which corresponding to wavelengths large compared to the effective size of the detector, 
a gravitational wave detector's antenna pattern has only monopole, quadrupole and 
sextupole components. Despite this limitation, it is in principle possible to locate a 
source with arbitrary precision so long as the signal-to-noise ratio is large and the 
motion of the detector with respect to the source is sufficiently fast and varied. 

Gravitational wave detectors get their directional information from the amplitude 
and frequency modulation that occur as the antenna pattern is swept across the sky. 
The frequency modulation is due to the Doppler shift caused by the relative motion 
of the source and the detector. For example, orbital motion about the Sun creates a 
periodic frequency shift with amplitude 5f/f ~ 10 _4 sin#, where 9 is the co-latitude in 
ecliptic coordinates. The amplitude modulation occurs as the antenna lobes are swept 
across the sky. For the LIGO detectors [p]] this occurs with a fundamental period of one 
sidereal day, while for the LISA detector the fundamental period is one sidereal year. 
When we are lucky enough to have two detectors, differences in the arrival time can be 
used to furnish further directional information. A discussion of the angular resolution of 
LIGO can be found in the work of Jaranowski & Kr61ak||. Discussions of the angular 
resolution of LISA can be found in the work of Peterseim et al. 0], Cutler || and Hellings 
& Moore©. 

In addition to individual bright sources of gravitational waves, there will also be a 
variety of unresolved gravitational wave backgrounds. It is likely that one component 
of the background will have a cosmological origin - the cosmic gravitational wave 
background (CGB). However, the dominant contribution is expected to be a confusion 
background formed by the superposition of many weak point sources, such as white 
dwarf binaries. The gravitational wave backgrounds will be approximately stochastic, 
making them difficult to distinguish from noise in the detectors. One way around this 
problem is to cross-correlate the output from two detectors. Since the noise in each 
detector is uncorrelated while the signal is correlated, the signal-to-noise improves as 
the square root of the observation time. Typically, this type of correlation has to be done 
over a finite frequency interval Af that is large compared to the Doppler modulation 
6f, so the frequency modulation is washed out. Thus, maps of the gravitational wave 
background have to be made using amplitude modulation alone. 

Here we describe how gravitational wave detectors can be used to create maps of the 
gravitational wave background. The pioneering work of Allen and Ottewill0 showed 
how anisotropies in the gravitational wave background will modulate the output from 
the cross-correlated LIGO detectors. A similar calculation was carried out by Giampieri 
and Polnorev|| in the low frequency limit (/ < 1 mHz) for the LISA detector. They 
showed how the signal would be modulated by the antenna's orbital motion. Neither 
group addressed the inverse problem - that is, how to create a map of the gravitational 
wave sky from the modulated output. We address the inverse problem and find that 
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the ground-based detectors that are currently under construction have the ability to 
produce maps with a resolution of I ~ 20 at / = 100 Hz, where I is the multipole 
number. This is comparable to the maps of the cosmic microwave background (CMB) 
produced by the COBE (COsmic Background Explorer) satellite. We find that a single 
LISA detector can produce a partial map with a resolution of i ~ 6 at / = 10 mHz, 
improving to £ ~ 25 at / = 100 mHz. More than one LISA detector needs to be flown 
to produce a complete map. For example, two LISA detectors with a relative orbital 
inclination of ~ 20 degrees could produce a complete map with a resolution of I ~ 8. 

A full implementation of the theoretical scheme developed here will be published 
elsewhere . It involves an end-to-end simulation of the data analysis, starting with 
a synthetic gravitational wave sky that produces a modulated signal from a simulated 
gravitational wave detector. The detector model incorporates a basic treatment of the 
anticipated noise. Finally, the simulated signal is deconvolved to produce a sky map 
that is compared to the original. 

2. The general convolution problem 

Suppose the gravitational wave background has a luminosity distribution given by 
£(9,<p,f). Here (9,4>) are coordinates on the sky in the CMB rest frame and / is 
the gravitational wave frequency. The response of a gravitational wave detector can 
be characterized by a detector response function T{9, 0, /, t). The response varies with 
time as the detector moves relative to the sky. The detector output C(t, f) is given by 
the all-sky integral 

c(tj) = I c(e,<i>j)r(e,(i>j,t) smodefy. (1) 

J s 2 

Given the time series C(t, f) we need to find the luminosity distribution C(9, 0, f). 

It is natural to decompose the luminosity distribution and the detector response 
function into spherical harmonics in their respective rest frames. Writing e a for basis 
vectors in the sky frame and e s for the detector frame we have 

e a = ts?Jt)ep A aP {t) = e a ■ ep , (2) 

where A a «(t) are the components of the coordinate transformation relating the two 
frames. In practice the gravitational fields will be weak and the velocities slow, so the 
coordinate transformation is well approximated by a spatial rotation. The coordinates 
(9, 0) in the sky frame are related to the coordinates (9, 0) in the detector frame by 

cos 9 = Q ■ e s = sin 9 cos A xZ + sin 9 sin A y2 + cos 9 A zS , (3) 

and 

sin 9 e 1 ^ = Q ■ (e s + ity) 

= sin 9 cos 4>{A xS . + iA X y) + sin 9 sin 4>(A yS + iA V y) + cos 9(A zS + iA Z y) (4) 
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We have written the transformations in this form since these are the combinations that 
appear in the spherical harmonics: 



2£ + l(£-m)\, , d 



(-sm6e^) m - — PAcose). (5) 



' \| 4tt (£ + m)\ V J d(cos9y 

Decomposing C(9, <p, /) and F(9, <j), f, t) in the sky frame we have 
£(e,(f>J) = Y,Pern{f)YUQ,<f>), and T(9, <f>, f,t) = ^a, m (/,t)y, m (^ <t>) , (6) 

im £m 

while in the detector's restframe we have 

F(6,t,f)=J2atm(f)Y im (eA). (7) 

The multipole moments in the two frames are related by 

i 

a >im(f>t) = X! ®tn(f)^imn(t) j (8) 

where 



A^ mn (t) = i y/ m (#, 0)y/„(e, 0) dn . (9) 

Using equations @, @ and @ it is a simple matter to calculate the transformation 
coefficients \£ mn (t). From the properties of the spherical harmonics it follows that 

K.mn= {- l ) m+n ^~m-n- (10) 

Thus, there we only need to calculate (21 + 1)(£ + 1) components at order I. Out to 
order I = 1 we have 



Ai_ 10 (t) = — (Kz(t) + iAy Z (t)) 



Aioi(t) = — j= (A zx (t) + iA zy (t)) 



Aooo(^) — 1 

1 

Ai_n(t) = - (-A xx (t) + A yy (t) - iA xy {t) - iA yx (t)) 
Aioo(t) = A zz (t) 

X no (t) = -j= (-A xz (t) + iA yz (t)) 
1 

Am(t) = - (A xx (t) + A yy (t) + iA xy (t) - iAy X (t)) . (11) 
Higher orders in £ involve higher powers of Aij(t). For example, 

A210W = \Kz{t) {-A xz {t) + zA^(t)) . (12) 

Using the decompositions defined in @, we can express the detector output in terms of 
the multipole moments pe m (f) and a^ m (t): 

C(t, f) = E(-!) m Ptm(f) f) ■ (13) 
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In most instances the detector scan pattern will be periodic in time, so it is natural to 
Fourier transform all time dependent quantities: 
1 r 27r 

A k = — dae- ika A(a), (14) 

Z7T JO 

where a = 2nt/T , and T is the period of the detector sweep. We can write 

C k(f) = Pim(f) IklmU) ( 15 ) 

where 

i 

lUmtf) = (-l) m a k£ . m (f) = (-I) 7 " E M/) ^ki-mn • (16) 

n=-l 

It is a simple matter to calculate the jk£m if we know the detector response function 
T and the scan pattern Ajj(t). The de-convolution problem comes down to solving the 
linear system of equations 

C k {f)=Pi{f)lki{f) where i = £ 2 + t + m + 1 . (17) 

The 7fcj(/) are known, the C k {f) are measured experimentally, and thepj(/) are what we 
want to find. Typically the system of equations will be both under and over constrained, 
but a best fit solution can be found by way of a singular value decomposition. 

3. Map making with LIGO 

We begin by studying how LIGO can be used to map the gravitational wave sky as 
a prelude to studying the space-based LISA detector. Allen and Ottewill|F| showed 
how the cross-correlated response of the Hanford and Livingston LIGO detectors is 
modulated as the rotation of the Earth sweeps the antenna pattern across the sky. 
However, they stopped short of solving the deconvolution problem to make a map of 
the gravitational wave sky. 

In a sky-fixed reference frame the response of a single LIGO detector to a plane 
gravitational wave h(t,x, /) with frequency /, propagating in the Q direction is given 
by 

s(t)=D(*):h(t,x), (18) 

where 

D(t) = - (u(t) ® u(t) - v(t) ® v(t)) (19) 

is the detector tensor and u(t) and v(t) are unit vectors in the direction of the 
interferometer arms. The response to a general gravitational wave h(t, x) can be derived 
from fll8|) by first decomposing the wave into a collection of plane waves: 

/°° r ^ ^ ^ ^ 

df dtt h A (f,tt)e- 2nift e 27TifQ -* /c e A (tt) . (20) 
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Here e + and e x are a set of polarization tensors. The signal we wish to analyze is formed 
by cross-correlating the output of the two LIGO detectors over an integration time T, 
centered at time t: 

r t+T/2 rt+T/2 

S(t) = alt' dt" Sl (t')s 2 (t")Q(t' -t"). (21) 

Jt-T/2 Jt-T/2 

Here Si(t) and s 2 (t) are the outputs of the two detectors and Q is some filter function. 
A stochastic background of gravitational waves can be characterized by the statistical 
character of the Fourier amplitudes ft,^/, Q). Following Allen and Ottewill we assume 
that the background can be approximated as a stationary Gaussian random distribution 
characterized by the expectation values 

(M/,«)) = o 

(h* A (fMAf'M = l -8{f - f) 5 -^p-5 AA , S h (f)P(Q) , (22) 

where Sh(f) is the power spectral density and P(£l) describes the angular distribution. 
From these relations it follows that the expectation value of S(t) is given by 

C(t) = (S(t)} = I r df S h (fh(f, t)Q(f) , (23) 

O J—OO 

where 

7 (/,t) = ^-/dnP(n) E F^t)F^t)e- 2 ^^' c . (24) 

87r J A=+,x 

Here Ax(t) is a vector connecting the Hanford and Livingston sites and 

F l A ({\,t) = V t (t):e A . (25) 

Using S h (f) = S h (-f), Q(f) = Q*(-f) and 7 (/,i) = >f{-f,t), we can write the 
signal as 

2T r°° / \ 

C{t) = — J o df S h (f) ( lR (f, t)Q R (f) - 7 /(/, t)Q!(f)) , , (26) 

where 7jj(/, t) denotes the real part of "f(f,t) etc.. In contrast to Allen and Ottewill, 
who considered the broad-band response, we are interested in making narrow band 
measurements over a small frequency interval of width Af centered at frequency /. 
This allows us to make maps of the gravitational wave sky at particular frequencies. To 
this end we use the top-hat filter 

QR + iqi f-¥<f'<f+%- 

QR-iqi -/-¥</'<-/ + ¥ (27) 

otherwise 
so that 



WO 



= 2 Sm(7rtA/) ( qR cos(2tt ft) + qj S m(2n ft)) . (28) 

Tit 

In the time domain the filter decays rapidly for times t > 1/A/. In choosing Af we 
have to ensure that 1/A/ T, otherwise the angular variation will be smeared out. 
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So long as Af is small compared to the range over which Sh(f) and j(f,t) vary, the 
narrow band response is given by 

C(t, f) = S h (f) (q RlR (f, t) - q ni (f, t)) . (29) 

This equation is of the form 

c(t,f) = <f ae,4>,f)r{e,4>,f,t) sinM?#, (30) 

where the luminosity equals 

c(e,<i>,f) = s h (f)P(fi), (31) 

and the detector response function is given by 

T% 0, /, t) = £ F^Cl, t)F 2 A (n, t) cos (2nfn • Ax(i)/c) 



+ IMSL V F^(h, t)Ff(n, t) sin (2vr/fi • Ax(t)/c) (32) 

The choice of integration period T results from a trade-off between maximizing the 
signal (T large) and maximizing the angular resolution (T small). Setting a maximum 
angular resolution of /c max fixes the integration period to be (Nyquist's theorem) 

T = ^. (33) 

It turns out that there is little to be gained by setting k max much larger than 50, so we 
will take /c max = 50 and set T = 861.64 s. Similar considerations apply to the choice of 
Af, which for LIGO we will take to be in the range Af = 1 — > 10 Hz. A plot of the 
antenna pattern for the cross-correlated pair of LIGO detectors is shown in Figure 1 for 
a range of frequencies. As expected, the antenna pattern has more structure at high 
frequencies. 

Now that we have the LIGO signal expressed in terms of the general convolution 
problem of the last section, it is a simple matter to carry through the analysis. We 
begin by setting a = 1-jitjT^ where Td = 86164 s is one sidereal day. Choosing Sky- fixed 
and Earth-fixed coordinate systems with their z axes aligned with the Earth's rotation 
axis, we find the two frames are related: 

e s = cos a e x + sin a e y 
Cy = — sin a e x + cos a e y 

e- z = e z , (34) 
where the bar denotes the Earth-fixed frame. From these relations it follows that 

^kimn = $mn $km) (35) 

and 

Ikim(f) = S km (-l) m ae- m {f) = Skmlemif) ■ (36) 



A=+,x 
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Figure 1. The antenna pattern for the cross-correlated LIGO detectors in the Earth- 
fixed frame, J 7 (9, <j>, /), with gi = and / = 25, 50, 100 and 200 Hz. 



Here jt m {f) correspond to the quantities calculated by Allen and Ottewill. However, 
they evaluated j£ m (f) in an Earth-fixed frame with the z axis parallel to the vector 
connecting the two detectors, Ax. A static rotation has to be applied to their results 
before we can use them. Once this is done we can attack the inversion problem 



C k(f) = J^PimifhMmif) = E PtkU)ltk(f) 
Im £=\k\ 

where the final line was obtained from the identity 

oo £ oo oo 

E = E E = E E • 



(37) 



(3f 



(in 



=0m=-< m=-oo £=\m\ 



At first sight it appears that we have no hope of solving the inversion problem since 
it involves 2(2fc max + 1) equations^] and an infinite number of unknowns. However, the 
ltm{f) decay rapidly for large £, effectively reducing the number of unknowns. The 
root-mean-square amplitude of the 7£ m (/)'s, 



Mf) = \JCt(f) , (39) 

provides a good estimate of how each order in I contributes to the sum in (0). Here 
Ci(f) is the usual angular power spectrum: 

1 



CtU) 



2£ + 



T E l7fon(/)| 5 

-L ™ a 



(40) 



% The equation count follows from — fc max < k < fc max and the freedom to alternately set qa and qi 
equal to zero. 
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Figure 2. The root-mean-square amplitudes of Jimif) for a range of frequencies. 



Plots of o~e(f) for / = 25, 50, 100 and 200 Hz are shown in Figure 2. In each case, the 
low multipoles dominate the response. Thus, the sum ( j37|) can be approximated by 



C k (f) 



(41) 



Plk{f)llk{f) , 
£=\k\ 

where £ max ~ 5 for / = 25 Hz, rising to ^ max ~ 18 for / = 200 Hz. The number of 
equations is restricted to those with \k\ < £ max - If we write C k = qnAu + qiB k and use 
two sets of filters, one with q R = 1, qj = 0, and the other with q R = and qj = 1, then 
(PD can be written as a set of 4£ max equations for (£ max + 1) 2 unknowns. As an example, 
consider the inversion problem at 25 Hz, where the system of equations is (for k > 0): 



A 


= P00700 + P20720 + P40740 


B 


= PlOllO + P30730 + P50750 


A x 


= £>2l721 + P41741 


B x 


= PXX1XX + P31731 + P51751 


A 2 


= P22722 + P42742 


B 2 


= P32732 + P52752 


A 3 


= P43743 


B 3 


= P33733 + P53753 


A 4 


= P44744 


B 4 


= P54754 






B 5 


= P55755 • 



(42) 
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The above system of 20 equations involves 36 p^ m 's, making it impossible to solve for each 
moment individually. The way around this impasse is to build another interferometer 
at a different location and to cross-correlate its output with the two LIGO detectors, 
in effect tripling the number of equations. Currently there are several interferometers 
being built in addition to the LIGO pair, including the Italian-French VIRGO detector 
near Pisa, the Anglo-German GEO detector near Hanover, and the Japanese TAMA 
detector near Tokyo. Cross-correlating the signals from all of these detectors will give 
roughly 40£ max equations for the (4na X + l) 2 unknowns, making it possible to create maps 
of the gravitational wave background with a resolution as high as £ max ~ 38. Indeed, 
the limitation on the resolution will not come from the system being under determined, 
but from the lack of antenna sensitivity at high I. 



4. Map making with LISA 

The procedure for making a map of the gravitational wave sky with the LISA detector 
is very similar to the procedure we used for the LIGO detectors. Currently there are 
no plans to fly two sets of LISA spacecraft, so we will probably have to make do 
with the self-correlated signal from a single detector. At first sight it would appear 
that correlated noise in the self-correlated signal will prevent us from being able to 
distinguishing detector noise from a stochastic gravitational wave signal. However, this 
turns out not to be the case. Firstly, it is possible to discriminate between detector noise 
and a stochastic signal by using the Sagnac signal that is formed by comparing the phase 
of a signal that is sent clockwise around the LISA triangle with a signal that is sent 
counter-clockwise pTI|, |T^ |. The Sagnac signal is very insensitive to gravitational waves, 
making it the perfect tool for monitoring instrument noise. Secondly, the signal will vary 
periodically as the detector sweeps across the sky whereas the noise will not. Thus, by 
making several complete sweeps, it is possible to build up the signal to noise ratio in 
the measurements of Ck for all k ^ 0. We study the map making capabilities of a single 
LISA interferometer and a pair of optimally cross-correlated LISA interferometers. The 
optimal cross-correlation is achieved by placing six LISA spacecraft in a circle, with two 
sets of three spacecraft forming independent interferometers rotated by 90 degrees in 
the plane of the circle . 



Proceeding as we did in the last section, the cross-correlated signal is given by 



C(t, f)=i C(9, 0, f)F(9, 0, /, t) sin 9dBd<j> , (43) 
Js 2 

where the luminosity is again given by 

C(9,cj ) ,f) = S h (f)P(n), (44) 
and the detector response function has the form 



Air 



27ri/n-Ax(t)/c 
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TAfqj 



Air 



£ F 1 A (n,f,t)Ff(n,f,t)e 2 ^ Q - A ^ c . (45) 

A=+,x 

While the basic expressions are similar for LISA and LIGO, the LISA antenna patterns 
are considerably more complicated: 

F l A (nj,t) = D l (nj,t):e A ({l) 

where 

1 



D(nj,t) 



(u(t) ® u(t)) r(u • n, f) - (v(t) ® v(t)) r(v • n, f) 



(46) 
(47) 



is the detector tensor and 



T(a-n,/) 



sine 



2/* 



1 - a • fi) ] exp ( -i^(3 + a • Q) 



f 



f 



+ sine ^(1 + a • fi) exp -irV(l + a • Q) 



2/ 



2/* 



(48) 



is the transfer function[|i"4]j. The transfer frequency, /* = c/ (2ttL), corresponds to a wave 
that just fits inside the interferometer. The LISA interferometer will have arms of length 
L = 5 x 10 6 km. The antenna pattern sweeps over the sky as the LISA constellation 
orbits about the Sun with period T = 1 sidereal year. Setting k mSuX = 50 fixes the 
integration period to be T = 315581 s (roughly three and a half days). 




Figure 3. The antenna pattern for the self-correlated LISA detector in the detector 
rest frame for a range of frequencies. 



Choosing Sky-fixed and LISA-fixed coordinate systems with their z axes normal to 
the ecliptic gives two frames of reference related by 

X X 
e s = —/=(2 — cos 2 a — cos a sina)^ H — ^(1 + cos 2 a — cos a sin a)e y 
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>/3, , c . w 
cos a + smaje 2 



e 5 = — t=(— 2 + cos 2 a — cos a sina)e x H — p(l + cos 2 a + cos a sin a)e y 



H — ^(cosa — smaje 2 



\/3 ^ \/3 ^ U 
e 5 = — cos ae^ + — sin ae^ + -e 2 , (49) 

where the bar denotes the LISA-fixed frame. The complicated scan pattern means that 
we have to calculate each rotation coefficient Xkimn individually. Since the coordinate 
transformation Ay (at) is second order in e ia , the index k will run from — 2£ to 2£. To 
order 1 = 1 the rotation coefficients are: 

Aoooo = 1 ^01-1-1 = (1 + i) ^11-10 — 

A21-11 = -jf- (1 — i) A_iio-i = — (1 + i) ^oioo = \ 

A1101 — (1 — *) A-211-1 — (1 + *) A-1110 — — ^ 

A0111 = (1 — i) A02-2-2 — f| A12-2-1 — (1 + i) 

^22-20 — 80 A 32 _21 — {1 — I) M2-22 ~ ~Jq 

A-12-1-2 = — Ai2_10 = A22-11 = (1 — i) (50) 

\ _ VSi \ — 3VT50i \ _ 3 /1 , -\ 

^32-12 — g— ^-220-2 — — 80 — A -120-l — ~g { l + 

1 \ 3/i -\ \ _ 3VT501 



^0200 — — 8 ^1201 — s (1 ~~ *) ^ 



2202 



^-321-2 — g - A -221-l — — I 1 + % ) ^-1210 — q~ 

Al212 = — A_ 422 _2 = A_ 3 22-l = — -JjP (1 + %) 

^-2220 — 80 A -1221 — — Jq U — % ) A 0222 — _ Jg 

Then antenna patterns for a self-correlated LISA detector and an optimally cross- 
correlated pair of LISA interferometers are shown in Figures 3, 4 and 5 for a range of 
frequencies in the rest frame of the detectors. These patterns can be turned into a^ m 's 
using the HEALPIX software package |I0[|. In combination with the analytically derived 



rotation coefficients, Xkimn, the a^ m 's yield the 7&| m 's that appear in the deconvolution 
problem (0). In the low frequency limit, / <C /*, it is possible to derive the antenna 
harmonics ae m analytically. For the self-correlated LISA detector the non-vanishing 
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harmonics at zero frequency are: 
a o = 1 
^20 - 

Cl4_4 



13 



2s/l 
7 



JO , -V210 
84 ~ 252 



(51) 



«40 



a 44 



12 



70 



V210 
252 



while for the optimally cross-correlated pair of LISA detectors 
a o = — 1 

a 20 — ~ 



n - V 70 
a 4 _4 — g^- 

n — ^ 
fl 44 — • 



(52) 
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Figure 4. The antenna pattern for two optimally cross-correlated LISA detectors in 
the detector rest frame for a range of frequencies. Here qn = 1 and qi = 0. 



Combined with the rotation coefficients Xke-mn, these give the 7jw m 's that we need 
to solve the deconvolution problem in the low frequency limit. For example, the self- 
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Figure 5. The antenna pattern for two optimally cross-correlated LISA detectors in 
the detector rest frame for a range of frequencies. Here qn = and qi = 1. 
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(53) 



Thus, at very low frequencies, the inversion problem only involves \k\ < 8 and 
multipoles with i = 0, 2 and 4, which gives us 17 equations for 15 unknowns. Writing 
out the convolution problem explicitly reveals that the system is simultaneously under 
and over constrained. In other words, some individual multipoles will be very well 
determined while others can only be given in combination with other multipoles. Once 
again the degeneracy can be split by building additional detectors that make different 
antenna sweeps. Alternatively, the LISA constellation could be re-position after several 
years of data collection into a new orbit that is inclined with respect to the ecliptic. 

At higher frequencies the LISA detectors are sensitive to more multipoles. Plots 
of the root-mean-square amplitude, ce(f), of the antenna multipoles ae m are shown in 
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Figure 6. The root-mean-square amplitudes of ae m (f) for the optimally cross- 
correlated LISA detectors at a range of frequencies. 



Figure 6 for the optimally cross-correlated LISA detectors. By either re-positioning the 
LISA constellation or by flying multiple constellations at different inclinations, it will 
be possible to produce maps with a resolution approaching £ = 25, which corresponds 
to an angular resolution of seven degrees. 
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